Method and system for measuring in a dynamic sequence of medical images 
Description 

Technical field 

The present invention relates to a method and a system for measuring in a 
dynamic sequence of medical images of a moving body part. 

Background of the invention 

Non-invasive methods for measuring functions and dysfunctions of body 
organs typically involve scanning with an X-ray device to obtain projection 
data for generation of a series of images. When measurements are to be 
performed of a moving body part, images of several events of the movement 
have to be created. Images are then generated with a certain frequency, 
whereby a sequence of consecutive Images is formed. The sequence depicts 
the same anatomical area over a period of time during which some dynamic 
event takes place and is thus a dynamic sequence. An example of such a 
dynamic sequence is an Image sequence of the heart over one or more 
heartbeat cycles. 

In order to be able to perform measurements In the images, different kinds 
of measurement tools for application to the images are needed. A 
fundamental tool in image measurements Is the distance measurement tool. 
Another important tool is the angle measurement tool for measuring the 
angle between two lines superimposed on an image. 

Common for many tools available today Is that they measure static data In 
static images. When a distance measurement is to be done in a static Image, 
a user Identifies two anatomical features in the image and applies a distance 
measurement tool to determine the distance between them. 

There are tools and methods for measuring in dynamic sequences, but they 
consist of two or more Instances of an intrinsically static tool that is more or 
less manually applied to several of the images in a dynamic sequence. The 
result from these different instances may then be used for comparisons or to 
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compute cruder approximations to dynamic entities by using ratios or 
differentiations. 

For example, a commonly used method today for measuring In a dynamic 
sequence Is a method of distance measurement applied to an image 
sequence from coronary angiography involving the following steps: manual 
identification of two images In the sequence, showing, for example, the end- 
diastole and end-systole phase in the same heart cycle; manual location of 
the most apical part of the coronary artery branches in both images; manual 
location of the lower contour of the left coronary ostium in both images; 
manual location of two points on the horizontal part of the circumflex artery, 
one proximal and one distal, in both images; measurement of the distance 
from the point from the second step to the point from the third step in one 
image and comparison with the distance between the same points in the 
other image; and measurement of the distance from the point from the 
second step to the points from the fourth step and comparison with the 
distance between the same points in the other image. 

Performing the above mentioned steps only once gives very limited dynamic 
information since only two images in the sequence are used for the 
measurements. To obtain more dynamic information, a user has to repeat 
the steps for every image in the dynamic sequence or at least for many 
images in the dynamic sequence. Usually there are over 100 images in a 
sequence. Thus, using the above mentioned method on all images in an 
image sequence is rather time-consuming. 

Through EP 1 088 517 a method and an apparatus for motion-free cardiac 
Imaging are known. In EP 1 088 517 a fixed reference point and a dynamic 
point are used for selecting images without motion-Induced artefacts in a 
30 dynamic sequence of computerized tomographic images. A line is drawn from 

the reference point to the dynamic point in all of the images. The line length 
represents the distance between the reference point and the dynamic point 
and one image or Images are selected in which the line length remains 
constant relative to the previous image. 

35 

Image recognition software may be used either to Identify a pair of reference 
points for measurements of relative motion or, once reference points are first 
Identified, to identify corresponding reference points on other Images. 
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Distances between the automatically identified points may be selected by 
software, based upon selected criteria. However, in EP 1 088 517 
determination of other variables such as speed, acceleration and retardation 
are not described for measurement in a dynamic sequence of images of a 
5 moving body organ. Furthermore, only two points, one fixed reference point 

and one dynamic point are used. It is an advantage to be able to use more 
than two points as well as more than one dynamic point. For example, angle 
measurements and area measurements require at least three points and 
measurements of relative motion of two moving parts of a moving body part 
10 require at least two dynamic points. 

Summary of the invention 

Accordingly, it Is an object of preferred embodiments of the present invention 
15 to provide an improved method for measuring in a dynamic sequence of 

medical images of a moving body part. 

This object my be achieved by means of a method for measuring in a 
dynamic sequence of medical images of a moving body part comprising: 

20 defining at least one measurement point in the moving body part in one of 

said images; defining a reference point in one of said images to a point being 
fixed relative to the image geometry, said reference point being different 
from said at least one measurement point; automatically tracking the at least 
one measurement point in all of said images of the sequence; automatically 

25 indicating the reference point in all of said images of the sequence; 

automatically determining a length and a direction of a vector extending from 
the reference point to one of the at least one measurement points for each 
pair of reference point and one measurement point in all of said images of 
the sequence, and automatically determining a rate of change of said length 

30 and/or said direction of said vector(s) between selected images in said 

sequence of images. 

Another object of preferred embodiments of the present invention Is to 
provide an improved method for generating a dynamic sequence of medical 
35 images of a moving body part and measuring in said dynamic sequence. 

This object may be achieved by means of a method for generating a dynamic 
sequence of medical images of a moving body part and measuring in said 
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dynamic sequence, which method comprises the steps of: scanning a portion 
of a body of a patient including the moving body part for generating time 
resolved projection data; generating said images from said projection data; 
defining at least one measurement point in the moving body part in one of 
5 said Images; defining a reference point in one of said images to a point being 

fixed relative to the image geometry, said reference point being different 
from said at least one measurement point; automatically tracking the at least 
one measurement point in all of said images of the sequence; automatically 
indicating the reference point in all of said images of the sequence; 

10 automatically determining a length and a direction of a vector extending from 

the reference point to one of the at least one measurement points for each 
pair of reference point and one measurement point in all of said images of 
the sequence, and automatically determining a rate of change of said length 
and/or said direction of said vector(s) between selected images in said 

15 sequence of images. 

A further object of preferred embodiments of the present invention is to 
provide an improved system for measuring in a dynamic sequence of medical 
images of a moving body part. 

20 

This object may be achieved through a system for measuring in a dynamic 
sequence of medical images of a moving body part, said system having 
means for: defining at least one measurement point in the moving body part 
in one of said images; defining a reference point in one of said Images to a 

25 point being fixed relative to the image geometry, said reference point being 

different from said at least one measurement point; automatically tracking 
the at least one measurement point in all of said Images of the sequence; 
automatically indicating the reference point in all of said images of the 
sequence; automatically determining a length and a direction of a vector 

30 extending from the reference point to one of the at least one measurement 

points for each pair of reference point and one measurement point in all of 
said images of the sequence, and automatically determining a rate of change 
of said length and/or said direction of said vector(s) between selected images 
in said sequence of images. 

35 

Another object of preferred embodiments of the present invention is to 
provide an Improved system for generating a dynamic sequence of medical 
images of a moving body part and for measuring in said dynamic sequence. 
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This object may be achieved through a system for generating a dynamic 
sequence of medical images of a moving body part and for measuring in said 
dynamic sequence, said system having means for: scanning a portion of a 
5 body of a patient including the moving body part for generating time resolved 

projection data; generating said images from said projection data; defining at 
least one measurement point in the moving body part in one of said images; 
defining a reference point in one of said Images to a point being fixed relative 
to the image geometry, said reference point being different from said at least 

10 one measurement point; automatically tracking the at least one 

measurement point in all of said images of the sequence; automatically 
indicating the reference point In all of said images of the sequence; 
automatically determining a length and a direction of a vector extending from 
the reference point to one of the at least one measurement points for each 

15 pair of reference point and one measurement point in all of said images of 

the sequence, and automatically determining a rate of change of said length 
and/or said direction of said vector(s) between selected Images in said 
sequence of images. 

20 Advantages of the methods and the systems according to preferred 

embodiments of the present invention include possibility to automatically 
measure movement variables such as speed, acceleration/retardation and/or 
direction of movement of a point or points of the moving body part, to use 
more than two points for measurement, to use more than one dynamic point 

25 for measurement, to automatically measure angles and areas and to use 

one-dimensional, two-dimensional and/or time resolved two-dimensional 
search fields for the automatically tracking of the dynamic point{s). 

Still other objects and advantages of the present invention will become 
30 apparent from the following detailed description considered in conjunction 

with the accompanying drawings. It Is to be understood, however, that the 
drawings are designed solely for purposes of illustration and not as a 
definition of the limits of the invention, for which reference should be made 
to the appended claims. It should further be understood that the drawings 
35 are not necessarily drawn to scale and that, unless otherwise indicated, they 

are merely Intended to conceptually illustrate the methods and systems 
herein. 
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Brief description of the drawings 

In the following, the invention is described in greater detail applied by way of 
example to a heart. The description is made with reference to attached 
5 drawings, in which like reference characters denote similar elements and 

Fig. 1 shows a composite of images in a dynamic sequence of schematic 
images of parts of a heart with a fixed reference point set in an apical part of 
the heart and a dynamic point set in a moving part of the left ventricular 
10 wall; 

Figs. 2a-d show separate images of the dynamic sequence shown in fig. 1 
and are also views of distance measurement; 

15 Figs. 3a-d show a first embodiment of the invention including the images 

shown in figs. 2a-d; 

Figs. 4a-d show a second embodiment of the invention including the images 
shown in figs. 2a-d; 

20 

Figs, 5a-d show the images shown in figs. 2a-d, but with a second dynamic 
point set in the left ventricular wall; 

Figs. 6a-d are views of measurement of dynamic angles and dynamic areas 
25 In the images shown in figs. 5a-d; 

Figs. 7a-d show a fourth embodiment of the present invention including the 
images shown in figs. 5a-d; 

30 Figs. 8a-d show a fifth embodiment of the present invention Including the 

images shown in figs. 5a-d; 

Figs. 9a-d show a sixth embodiment of the present Invention including the 
images shown in figs. 5a-d; 

35 

Fig. 10 shows a schematic view of four points used for angiographic 
measurements. 
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Description of embodiments 

Figure 1 shows a composite of highly schematic images of some parts of a 
heart 1. The images shown in figure 1 belong to a dynamic sequence of 
5 Images of the heart 1, but are only a few images of a dynamic sequence 

covering a cardiac cycle of the heart 1. The sequence of images may be 
generated from projection data from time resolved two dimensional X-ray 
scanning of a portion of a body of a patient including the heart 1. The 
projection data may thus be time resolved two-dimensional data. In order to 
10 cover a complete cardiac cycle the images may be generated at an adequate 

frequency, for example 12,5 images per second. A contrast medium may be 
introduced before the scanning to visualize vessels properly on the Images. 

The images could be generated at any frequency that permits the analysis of 
15 the relevant moving body part. For example, the Images could be generated 

at a rate of from about 9 Images per second to about 24 Images per second. 
However, the rate could be greater than about 24 images per second, 
particularly for fast moving body parts. On the other hand, the rate could be 
less than 9 Images per second, particularly for slow moving body parts. 

20 

The heart 1 is a cyclically moving body part having a right atrium 2, a left 
atrium 3, a right ventricle 4 and a left ventricle 5. A single cardiac cycle 
consists of one diastolic phase (expansion phase) and one systolic phase 
(contraction phase) for the atria 2, 3 and one diastolic phase and one systolic 

25 phase for the ventricles 4, 5. The right ventricle 4 receives deoxygenated 

blood during its diastolic phase and pumps it into the pulmonary arteries 
during its systolic phase, while the left ventricle 5 receives oxygenated blood 
during Its diastolic phase and pumps it into the aorta during its systolic 
phase. The heart muscles are supplied with blood by the coronary arteries, 

30 which encircle the heart and thus follow the expansion and contraction of the 

heart 1 during a cardiac cycle. 

For different reasons it would be advantageous to be able to track a point of 
the heart 1 and measure for example movement variables of that point of the 
35 heart 1 relative a fixed point during the cardiac cycle. Movement variables, 

which might be interesting to measure may Include distance, speed, 
acceleration, retardation, direction of movement of a point of the heart as 
well as dynamic angles and areas. 
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Using the method according to a preferred embodiment of the invention, a 
reference point 6 and a measurement point 7 typically are defined by a user. 
The reference point 6 may be set to a point that is fixed in its position 
relative the Image geometry. That is, the reference point may be set to a 
point of a structure that is not essentially moved during a cardiac cycle. In 
the embodiment illustrated by the images shown in figure 1, the reference 
point 6 is set to an apical part of the heart 1, which Is not essentially moved 
during a cardiac cycle. The reference point 6 Is set in one image by a user 
and is then automatically Indicated, that is, marked, in the other images in 
the sequence through Image processing software. 

In the embodiment Illustrated by the images in figure 1, a point of the wall of 
the left ventricle 5 Is desired to be measured and then that point is defined 
as the measurement point 7. The measurement point 7 may be a dynamic 
point, since it is moved during the cardiac cycle. The dynamic point 7 may 
also be set by a user in one image, and the corresponding points in the other 
images may then be automatically tracked by using image processing 
software including one of a number of known algorithms for automatic 
tracking of anatomic parts. An example of one such algorithm is FMI-SPOMF 
(Fourier-Mellin Invariant Symmetric Phase-Only Matched Filtering), described 
for instance by Qin-sheng Chen, Michel Defrise and F. Deconick in IEEE 
Transactions on Pattern Analysis and Machine Intelligence, Vol. 16, No. 12, 
pp 1146-1168, Dec 1994. In order to improve the efficiency of the 
algorithm/algorithms used, the images may be pretreated to increase the 
contrasts between the imaged object and the background. The tracking may 
be performed using information about the imaged object and the darkest 
point may, for example, be tracked since it probably represents a vessel 
containing a lot of contrast medium. 

The wall of the left ventricle 5 expands during the diastolic phase and 
contracts during the systolic phase. Thus, during the diastolic phase of the 
left ventricle 5, the dynamic point 7 typically moves from a position denoted 
8 into expanded positions denoted 9, 10 and 11, The position 8 represents 
the end-systolic position of the point 7 and the position 11 represents the 
end-dlastolic position of the point 7. Accordingly, during the systolic phase of 
the left ventricle 5, the point 7 may move from the position 11 into position 
8. 
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Figures 2a-d show separate images of the images in the dynamic sequence 
shown in fig. 1 and are also views of distance measurement. The images in 
figures 2a-d are subsequent in time. Figure 2a shows an end-systolic view of 
5 the heart 1 and figure 2d an end-diastolic view of the heart 1, A length and a 

direction of a vector Vi extending from the reference point 6 to the first 
dynamic point 7 may be automatically determined in all of the images of the 
sequence using image processing software. A first distance du between the 
dynamic point 7 and the reference point 6 may also be automatically 
10 determined in each of the images in the sequence of images using the length 

of the vector Vi. The first distance du between the reference point 6 and the 
dynamic point 7 may also be automatically compared between different 
images and used to determine the change of the position, i.e. the movement, 
of the dynamic point 7 during the cardiac cycle. 

15 

Since the data used for reconstruction of the images in figure 1 may be time 
resolved two-dimensional data, it may also be possible to determine a rate of 
change of the first distance du, i.e. movement variables, between different 
images in the sequence of images. The rate of change may also be 
20 automatically determined using image processing software. 

Movement variables which may be determined may include for example 
speed, acceleration, retardation and direction of movement of the dynamic 
point 7 in different moments of the cardiac cycle, i.e. in different images of 

25 the sequence of images. Furthermore, the speed, acceleration and/or 

retardation with which the dynamic point 7 is moved between its different 
positions in two different images may also be determined. The movement 
variables may be measured in all of the below described embodiments of the 
present invention. When determining the rate of change of movement, the 

30 rate of change of any one or more of the variables may be determined. 

Figures 3a-d show a first embodiment of the invention including the images 
shown in figures 2a-d. In the first embodiment of the invention a one- 
dimensional search field may be used for the tracking of the dynamic point 7 
35 in the different images in the sequence of images. A vector 12, i.e. a line, 

representing the one-dimensional search field may be marked in one of the 
images and the dynamic point 7 may then be automatically tracked along 
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that line 12 in all of the Innages in the sequence of innages using image 
processing software as previously mentioned. 

Figures 4a-d show a second embodiment of the invention including the 
images shown in figures 2a-d. In the second embodiment of the invention a 
two-dimensional search field is used for the tracking of the dynamic point 7 
in the different images in the sequence of images. An area 13, i.e. a 
rectangle, representing the two-dimensional search field may be marked in 
one of the images and the dynamic point 7 may then be automatically 
tracked within that rectangle 13 in all of the images in the sequence of 
images using image processing software as previously mentioned. 

In a third embodiment (not shown) of the invention a time resolved two- 
dimensional search field may be used for the tracking of the dynamic point 7 
in the different images In the sequence of images. An area 13, I.e. a 
rectangle, representing the two-dimensional search field may be marked in 
one of the images in the sequence of images and the dynamic point 7 may 
be automatically tracked within the rectangle using image processing 
software as previously mentioned, but information from previous and 
subsequent images may be used as guide to which point of possible points in 
the rectangle that is the searched dynamic point 7. Expected positions of the 
dynamic point 7, which may be based on information from previous images, 
may also be used to choose which point in the rectangle that is the searched 
dynamic point 7. 

Figures 5a-d show the images shown in figures 2a-d, but with a second 
dynamic point 14. Using the method according to the invention, movement of 
two different parts of the moving body part may simultaneously be measured 
by defining a second dynamic point 14, which may be different from the first 
dynamic point 7. In the figures 5a-d, the second dynamic point 14 may be 
set to a point of the wall of the left ventricle 5 being different from the first 
dynamic point 7. The second dynamic point 14 may be moved from a 
position denoted 19 into expanded positions 20, 21 and 22. The position 19 
may represent the end-systolic position of the point 14 and the position 22 
may represent the end-diastolic position of the point 14. Accordingly, during 
the systolic phase of the left ventricle 5, the point 14 may be moved from 
position 22 into position 19. A length and a direction of the vector Vi 
extending from the reference point 6 to the first dynamic point 7 as well as a 
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vector V2 extending from the reference point 6 to the second dynamic point 
14 may then be automatically determined using image processing software. 

The first distance dn between the first dynamic point 7 and the reference 
point 6 as well as the first distance between the second dynamic point 14 
and the reference point 6 may then also be automatically determined using 
the length of the vector Vj and the vector Vz, respectively, and compared 
between different Images. Furthermore, the movement between the position 
in one image and the position in another Image of both the first dynamic 
point 7 and the second dynamic point 14 as well as movement variables like, 
for example, the above mentioned may also be determined. It is also 
possible to automatically determine the movement of the two dynamic points 
7, 14 relative each other through determining a second distance dzi between 
the two dynamic points 7, 14 and comparing the second distance d2i 
between different Images. The second distance dzi may be automatically 
determined using the lengths of the vectors Vi and V2 and using Image 
processing software. Furthermore, using the two dynamic points 7, 14 it is 
also possible to measure dynamic angles and areas. Dynamic angles and 
areas may be measured In all of the below described embodiments. 

Figures 6a-d are views of measurement of dynamic angles and dynamic 
areas in the images shown in figures 5a-d. A dynamic angle x may be 
automatically determined by using image processing software and by using 
the first distances dn, d^ and the second distance d2i in each of the images 
of the sequence. Other angles than the shown angle x may of course also be 
determined. A dynamic area a may automatically be determined by using 
Image processing software and by using the first distances dn, d^ and the 
second distance d2i in each of the images of the sequence. 

30 Figures 7a-d show a fourth embodiment of the invention including the images 

shown in figures 5a-d. In the fourth embodiment of the invention one- 
dlmenslonal search fields may be used for the tracking of the first and second 
dynamic points 7, 14, respectively. In the different Images in the sequence of 
images. A vector 12, i.e. a line, representing the one-dimensional search field 

35 may be marked for each dynamic point 7, 14 in one of the Images and the 

dynamic points 7, 14 may then be respectively tracked along the lines 12 In 
all of the Images in the sequence of images using Image processing software 
as previously mentioned. 
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Figures 8a-d show a fifth embodiment of the invention including the images 
shown in figures 5a-d. In the fifth embodiment of the invention two- 
dimensional search fields may be used for the tracking of the first and second 
5 dynamic points 7, 14, respectively, in the different images In the sequence of 

images. An area 13, i.e. a rectangle, representing the two-dimensional 
search field may be marked for each dynamic point 7, 14 in one of the 
images and the dynamic points 7, 14 may then be respectively tracked within 
those rectangles 13 in all of the images in the sequence of images using 
10 image processing software as previously mentioned. 

Figures 9a-d show a sixth embodiment of the invention including the images 
shown In figures 5a-d. In the sixth embodiment of the invention a one- 
dimensional search field may be used for the tracking of one of the first and 
the second dynamic points 7, 14 and a two-dimensional search field may be 
used for the tracking of the other of the first and the second dynamic points 
7, 14 in the different images in the sequence of images. A vector 12, i.e. a 
line representing the one-dimensional search field and an area 13, i.e. a 
rectangle, representing the two-dimensional search field may be marked in 
one of the images and the dynamic points 7, 14 may then be tracked along 
the line and within the rectangle, respectively. In all of the images In the 
sequence of images using Image processing software as previously 
mentioned. In the figures 9a-d a one-dimensional search field may be used 
for the tracking of the first dynamic point 7 and a two-dimensional search 
field may be used for the tracking of the second dynamic point 14, but 
instead a one-dimensional search field may be used for the tracking of the 
second dynamic point 14 and a two-dimensional search field may be used for 
the tracking of the first dynamic point 7. 

30 In a seventh embodiment (not shown) of the invention time resolved two- 

dimensional search fields may be used for the tracking of the dynamic points 
7, 14, respectively, in the different images in the sequence of images. An 
area 13, i.e. a rectangle, representing the two-dimensional search field may 
be marked for each dynamic point 7, 14 in one of the images in the sequence 

35 of images and the dynamic points 7, 14 may be respectively tracked within 

the rectangles using image processing software as previously mentioned, but 
information from previous and subsequent images is used as guide to which 
points of possible points in the rectangles that are the searched dynamic 
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points 7, 14. Expected positions of the dynamic points 7, 14 based on 
information from previous images may also be used to choose which points 
that are the searched dynamic points 7, 14. 

5 In an eighth embodiment (not shown) of the Invention a one-dlmenslonal 

search field may be used for the tracking of one of the first and the second 
dynamic points 7, 14 and a time resolved two-dimensional search field may 
be used for the tracking of the other of the first and the second dynamic 
point 7, 14 in the different Images in the sequence of Images. A vector 12, 

10 I.e. a line representing the one-dimensional search field and an area 13, i.e. 

a rectangle, representing the two-dimensional search field may be marked in 
one of the Images and the dynamic points 7, 14 may then be respectively 
tracked along the line 12 and within the rectangle 13 in all of the images In 
the sequence of Images using Image processing software as previously 

15 mentioned. However, information from previous and subsequent images may 

be used as guide to which point of possible points in the rectangle 13 that 
may be the searched dynamic point 7, 14. Expected positions of the dynamic 
point 7, 14 based on information from previous Images may also be used to 
choose which point in the rectangle that Is the searched dynamic point 7, 14. 

20 

In a ninth embodiment (not shown) of the invention a two-dimensional 
search field may be used for the tracking of one of the first and the second 
dynamic points 7, 14 and a time resolved two-dimensional search field may 
be used for the tracking of the other of the first and the second dynamic 

25 points 7, 14 in the different images In the sequence of images. An area 13, 

i.e. a rectangle, representing the two-dimensional search field may be 
marked for each dynamic point 7, 14 In one of the images and the dynamic 
points 7, 14 may then be respectively tracked within that rectangles 13 In all 
of the Images In the sequence of Images using image processing software as 

30 previously mentioned. However, information from previous and subsequent 

images may be used as guide to which point of possible points in the 
rectangle 13 representing the time resolved two-dimensional search field that 
is the searched dynamic point 7, 14. Expected positions of the dynamic point 
7, 14 based on information from previous images may also be used to choose 

35 which point In the rectangle 13 representing the time resolved two- 

dimensional search field that may be the searched dynamic point 7, 14. 
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The method according to the invention may be expanded to cover more than 
two dynamic points as well as more than one reference point. Any 
combinations of one-dimensional, two-dimensional and time resolved two- 
dimensional search fields for the tracking of the different points may then be 
possible to use. 

By way of example, an embodiment of the method according to the invention 
is described involving coronary angiography. Projection data of the coronary 
arteries are generated through coronary angiography, i.e. X-ray examination 
after injection of an opaque dye. The Images are generated during at least 
one cardiac cycle at a frequency of 12.5 images per second. The right 
anterior oblique 30** projection has shown good results and is preferably 
used. A sequence of images may then be generated through reconstruction 
from the projection data. Four points 15, 16, 17, 18 may then be located and 
marked In one Image in the sequence of images; a point 15 representing the 
most apical point of the coronary artery branches, a point 16 representing 
the lower contour of the left coronary ostium, a proximal point 17 and a 
distal point 18 on the horizontal part of the circumflex artery. The 
corresponding four points 15, 16, 17, 18 may then be automatically located 
in each of the images In the image sequence using image processing 
software. Figure 10 shows a schematic view of the four points 15, 16, 17, 18 
used for angiographic measurements. 

The epicardial part of the apex is almost stationary during the heart cycle. 
Thus, the most apical point 15 of the coronary artery branches is almost 
stationary during the heart cycle and may therefore be used as a reference 
point. The other three previously mentioned points 16, 17, 18 may be 
dynamic points representing a first, a second and a third dynamic point. A 
one-dimensional, a two-dimensional or a time resolved two-dimensional 
search field may be used for the tracking of each of the corresponding 
dynamic points 16, 17, 18 in the other images in the sequence of images. 
The first distance du from the most apical point 15 of the coronary artery 
branches to the point 16 representing the lower contour of the left coronary 
ostium may then be automatically determined in each image showing the 
amplitude of the left coronary ostium motion. The first distance d^ between 
the most apical point 15 of the coronary artery branches and the proximal 
point 17 on the horizontal part of the circumflex artery as well as the first 
distance di3 between the most apical point 15 of the coronary artery 
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branches and the distal point 18 on the horizontal part of the circumflex 
artery may also be automatically determined in each image showing the 
motion of the proximal point 17 and the distal point 18, respectively. 
Furthermore, the second distance d223 between the proximal point 17 and the 
distal point 18 of the circumflex artery may be determined, which second 
distance d223 represents a portion of the circumflex artery extending along 
the atrioventricular groove. The motion of that portion of the circumflex 
artery roughly follows the most basal part of the left ventricular wall, which 
means that change of the distance between the proximal and the distal point 
of the circumflex artery represents the shortening and lengthening of the left 
ventricle 5. 

Although cardiac imaging embodiments are described in detail, the invention 
has more general applicability and may be applied to another moving body 
part, such as for example the lungs. Furthermore, the method according to 
the invention may also be expanded to work on time resolved three- 
dimensional angiographies as well as simultaneous measuring of dynamic 
points in different body parts. 

Thus, while there have been shown and described and pointed out 
fundamental novel features of the invention as applied to embodiments 
thereof, it will be understood that various omissions and substitutions and 
changes in details of the methods described, and in their operation, may be 
made by those skilled in the art without departing from the spirit of the 
invention. For example, it is expressly intended that all combinations of those 
method steps and/or system elements which perform substantially the same 
function in substantially the same way to achieve the same results are within 
the scope of the invention. I^oreover, it should be recognized that method 
steps and/or system elements shown and/or described in connection with 
any disclosed form or embodiment of the invention may be incorporated in 
any other disclosed or described or suggested form or embodiment as a 
general matter of design choice. 



